Many-body physics of intersubband polaritons 
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Intersubband polaritons are light-matter excitations originating from the strong coupling between 
an intersubband quantum well electronic transition and a microcavity photon mode. In this paper 
we study how the Coulomb electron-electron interaction and the Pauli saturation of the electronic 
transitions affect the physics of intersubband polaritons. We develop a microscopic many-body the- 
ory for the physics of such composite bosonic excitations in a microcavity-embedded two-dimensional 
electron gas. As a first application, we calculate the modification of the depolarization shifts and 
the efficiency of intersubband polariton-polariton scattering processes. 



I. INTRODUCTION 

Intersubband polaritons are elementary excitations of 
microcavities embedding doped semiconductor quantum 
wells. In particular, they are the normal modes due to 
the strong coupling between a cavity photon mode and 
the transition between two conduction subbands in quan- 
tum wells containing a two-dimensional electron gas. In- 
tersubband polaritons were experimentally demonstrated 
for the first time in 2003^ and their physics is the object 
of many interesting investigations, both fundamental^— 
and applied^ 3 - - — . 

So far, theoretical investigations have been based 
mostly on a bosonized Hamiltonian approach 2 . In the 
dilute regime, i.e., when the excitation density is much 
smaller than the electronic density in the doped quan- 
tum wells, intersubband polaritons behave as (compos- 
ite) bosons. At higher excitation densities anyway, Pauli 
saturation effects due to electron fermionic statistics and 
deviations from the bosonic behavior are expected 2 ^. 

An important aspect of the physics of such excitations 
is the role of the Coulomb interaction. In the case of bare 
intersubband transitions, Coulomb interaction is known 
not to modify dramatically the polaritonic lineshape^ 3 -, 
being responsible mainly for a simple shift of the exci- 
tation frequency, the so-called depolarization shift. This 
is why, in many relevant situations, the presence of the 
Coulomb interaction can be in first approximation ne- 
glected, by using a renormalized (measured) value for 
the frequency of the intersubband transition instead of 
the bare one, calculated from the single-electron subband 
structure. 

Very recently, the role of Coulomb interaction in inter- 
subband polariton systems has started to attract interest. 
In Ref. [2J] the part of the Coulomb interaction respon- 
sible for the depolarization shift has been studied in the 
bosonic excitation formalism, allowing to calculate the 
renormalization from first principles. Other works have 
studied the same problem using different gauge s 25 ! 26 and 
a Green functions approach 2 ^. As a matter of fact, up to 
now, only depolarization or static effects have been taken 
into account. Hence, the general problem of the many- 



body physics of intersubband polaritons is largely unex- 
plored, in particular the physics of polariton-polariton 
interactions and polariton-polariton scattering. 

In this paper we derive a microscopic many-body the- 
ory describing the intersubband polariton-polariton scat- 
tering due both to the Coulomb interaction and to the 
non-perfect bosonicity of the intersubband excitations. 
As a first application of the present theoretical formalism, 
we calculate how the depolarization shifts are modified 
and, most importantly, we determine the efficiency of the 
intersubband polariton-polariton scattering processes. 

This article is organized as follows. In Sec. HH we 
present the physical system, introduce the fermionic 
electron-photon Hamiltonian including the different 
channels of the (screened) Coulomb interaction. More- 
over, we write the intersubband excitations in second 
quantized formalism. In Sec. IIII1 the composite boson 
commutator approach 28 to calculate many-body matrix 
elements is presented for the case of intersubband excita- 
tions. In Sec. IIV1 an effective bosonic Hamiltonian (up 
to the fourth order in the bosonic operators) is derived 
in such a way that the two-excitation properties are the 
same as those predicted by the microscopic composite 
boson approach. In Sec. IVl numerical applications of the 
present theory to the calculation of the interaction energy 
between intersubband excitations are presented. Sec. 
IVII is devoted to the calculation of polariton-polariton 
scattering processes. Finally, conclusions are drawn in 
Sec. I VIII Additional technical details are included in 
Appendix \Z\ E and 



II. THE SYSTEM AND ITS HAMILTONIAN 
DESCRIPTION 

A. Physical system 

We will consider a photonic, planar microcavity em- 
bedding nQw symmetric quantum wells, each doped with 
an areal electron density of n e i (see Fig. [T]for a schematic 
representation of such a system). As the quantum wells 
confine the electrons along the growth direction, the elec- 
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FIG. 1: Schematic representation of the system under inves- 
tigation. Panel (a): Multiple doped quantum wells (QWs) 
are embedded in a planar photonic microcavity. The image 
above shows also the fundamental TM photonic mode in the 
case of metallic mirrors. Due to intersubband selection rule, 
the intersubband transitions only couple to the z-component 
of the electric field. Panel (b): Energy profile of the quan- 
tum well along the z-direction and its two bound states (left). 
Electronic dispersion of the first two conduction subbands in 
each of the doped quantum wells, as a function of the in- 
plane electronic wavevector (right). The Fermi wavevector, 
kF, the Fermi Energy, Ef, and the intersubband gap, ftto\2, 
are highlighted. 



tronic bands will be split into multiple, almost parallel 
subbands. The optically active transitions we are inter- 
ested in will thus be the ones from the last occupied 
conduction subband to the first unoccupied one. The 
photonic microcavity width will be chosen to have one of 
the photonic modes almost resonant with the electronic 
transition between these two electronic conduction sub- 
bands (the first two for simplicity, even if more general 
situations are possible^). In order to simplify the no- 
tation, we will consider that the quantum wells are all 
identically coupled to the relevant microcavity field. 

The Hamiltonian describing the coupled electron- 
photon system can be written in the Coulomb gauge as 



H 



He + Hf 



(1) 



where He includes only photonic degrees of freedom, 



while Hp regroups all the terms presenting also electronic 
ones. Introducing the annihilation operators c Mi k for an 
electron in the subband /i with an in-plane wavevector 
k, and a q for a microcavity photon with an in-plane 
wavevector q, we can write the two Hamiltonians explic- 
itly as 



H c = 



H F 
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E^' 

q 

-^Kin 



/jA g (a q + aL q )(a_ q + a q ) 



H] 



Rabi 



H, 



Coul, 



(2) 



ftuj 



H 



Rabi — 



Hi 



Coul — 77 



E 

^={l,2},k 

E^ C 2,k- -q 
k,q 

1 

2 

/i,,/i/,i/,i/£{l,2} 



(3) 



c i,k + c i,k+q c 2,k)(a q + aL q ) 



E 

k,k',q 



v q V,k+q°l/,k'- 



t^V' ,k' C/i' ,ki 



are respectively the electronic kinetic energy, the light- 
matter coupling Hamiltonian and the electron-electron 
Coulomb interaction. Coefficients w Ci9 and are the 
frequencies of a photon with in-plane wavevector q and 
of an electron with in-plane wavevector k in the subband 
/i, Xq is the light-matter coupling frequency, A 9 the dia- 
magnetic interaction coefficient and the q = term has 
to be excluded from the sum in Hcoui in order to com- 
pensate for the homogeneous positive background of the 
jcllium models. 

In order to simplify the notation both spin and quan- 
tum well indexes have been omitted but one should keep 
in mind that there is an implicit sum over them. The 
explicit notation is given in Appendix [SJ We also neglect 
to explicitly write the photon polarization because, due 
to the intersubband selection rule, only the TM modes 
are coupled with matter. 

The Coulomb matrix elements V^ vv M are given by^I 



q 



q 



2Seoe r q 9 



(4) 



where S is the surface of the sample and ^ s ^ ne rea l 

^-component wavefunction of an electron in the subband 

Due to the symmetry of the quantum wells, and thus of 
the wavefunctions, a certain number of matrix elements 
in Eq. Q can be seen to be zero, in particular all matrix 
elements with an odd number of I and 2 indices 
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The other elements obey the relations 
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so that there are only four distinct matrix elements, 
V q 1212 ,V q 1221 ,V q lln and V 2222 . In Fig. H we show a graph- 
ical representation of these four qualitatively different 
processes and the wavevector dependency of their re- 
spective coefficients for an infinite potential well. From 
such a figure we can clearly see a qualitative difference 
in the long wavelength behavior between V q 1212 and the 
other three coefficients. As explained in greater detail 
in Ref. [Hj], they indeed refer to qualitatively different 
physical processes. While V q 1212 describes an intersub- 
band, coherent Coulomb scattering, that is responsible 
for the so-called depolarization shift, the others describe 
intrasubband scattering processes, in which electrons are 
scattered inside the same subband. 



B. Screened Coulomb interaction 

The dense two dimensional electron gas in the first 
electronic subband screens the Coulomb interaction. In 
order to take it into account we will thus replace the bare 
Coulomb interaction V q avu M defined by Eq. d3|, with its 

static RPA-screened version y-fv'S, FromRefs. [HH, 
these coefficients obey the Dyson equation 



v q 



(7) 



where II Q ^(q, u>) is the RPA polarization function for the 
a — > f3 transition. The screened interactions thus take 
the form 



V, 



V: 



V, 
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(8) 
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V 
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i-y (/ 1212 [n 12 (q,o) + n 21 (q,o)]' 



where v = {1,2}. Screening of intrasubband and inter- 
subband processes can thus be encoded respectively in 
the dielectric functions e(q) and ei2(q), where 



and 



e (q) = l-y g ull n 11 (q ) 0), (9) 



£12 (q) = 1 - K 1212 [n 12 (q,0) + n 21 (q,0)]. (10) 



As we will see in the rest of this Section, only the long- 
wavelength limit is relevant in the study of intersubband 
polaritons. We will thus always consider the q — > limit 
of the above experessions. In Fig. [2] the reader can find 
the screened and unscreened profiles of the different po- 
tentials as a function of the exchanged momentum nor- 
malized over the Fermi wavevevector q/kp. 



C. Intersubband excitations 

The theory in this paper will be developed using a zero 
temperature (T = 0) formalism. In order for this approx- 
imation to be valid, we need to have a vanishing ther- 
mal electronic population in the excited subband, that 
is hui2,o ~~ Ep 3> ksT, where Ep is the Fermi energy 
and ks the Boltzmann constant. Under such a condition 
the ground state of the uncoupled light-matter system is 
given by 



w - n 



j 

"l,k 



|o> 



(11) 



k|<fc F 



where kp is the Fermi wavevector and |0) is the vacuum 
state, with no electrons in both subbands and no photons 
in the microcavity. In the ground state all electronic 
states up to the Fermi energy are occupied and all the 
others are empty. 

If we neglect the incoherent intrasubband Coulomb 
processes (due to the coefficients V 2222 , V q 1221 and V q ln ) 
in Eq. the full Hamiltonian H can be solved, in the 
dilute regime, in terms of bosonic coherent excitations 
named intersubband polaritons 2 . 

In order to do this, we first notice that the microcav- 
ity photon mode is coupled through i?R a bi to only one 
particular linear superposition of electronic transitions 



i 



E* t 
^0,k C 2,k+q 



Cl,k, 



(12) 



where ^o,k = Q(kp — k) and is the Heaviside func- 
tion. The choice of the notation with the index will 
be explained in the following. In Eq. (jT2"j) . as well as in 
Eq. ([3]), the sum is implicit over the spin and quantum 
well indexes to simplify notations. This defines an inter- 
subband collective excitation, which we will name bright 
intersubband excitation. Hamiltonian i?R, a bi can thus be 
written exactly 



(13) 



where il q = ^/nQwNc\Xq is th e Rabi frequency. Eq. (fT3|) 
makes it clear that bright excitations concentrate all the 
oscillator strength of the electron gas. Other similar col- 
lective excitations can be constructed by an orthonormal- 
ization procedure 
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•QwAcl 



k 



' y i,k c 2,k+q c l.k> 



(14) 



where index i runs from 1 to UQw-^el — 1 an d the v coef- 
ficients have support over the Fermi sea and satisfy the 
orthonormality relation 
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FIG. 2: Panels in the first two lines show the wavevector dependency of the the four qualitatively different Coulomb processes. 
In each panel we plot both the bare potentials (dashed line) and the static-RPA screened ones (solid line). The first row is 
for a mid-infrared transition, with huji2 = 140meV and an electronic density in each quantum well n e i = 5 x 10 11 cm -2 . The 
second row is instead for a THz transition, with Hui2 = 15meV and n e i = 2 x 10 11 cm" 2 . In the third row we present schematic 
representation of the four relevant processes. 



However, none of these new collective excitations is cou- 
pled to the microcavity photon field. This is why we call 
them dark intersubband excitations. We will consider 
them in Sec. IIIII 

We then use the fact that the two electronic subbands 
are parallels, and we neglect the photonic exchanged mo- 
mentum when compared to an electronic one, that is 

w 2,fc = Vl,k + u 12 
W |i,|k+q| - W^fc. (16) 

We can thus rewrite the Hamiltonian 

H = lie I X^' A../'.- '! ( 17 ) 

+ ^g(60,-q + ^, q )(°q + °-q) 

q 

+^En m2 K + & o,- q ) (4,- q + ^o, q ), 
q 

where Hq is defined in Eq. ^ and the last line comes 
from the intersubband part of the Coulomb interaction 



only (we have omitted the constant ground state energy). 
If the &o,q operators are supposed to obey bosonic com- 
mutation rules 

[ 6 o,qX,q'] - S wSi,j, (18) 

then the Hamiltonian in Eq. f| 1 T[) is quadratic 
and can thus be solved by an Hopfield-Bogoliubov 
transformation 3 ^. 

The different approximations that led to Eq. (fT7|) 
(parabolicity, small exchanged momentum, perfect 
bosonicity, absence of intrasubband scattering) have been 
thorougly tested in multiple experiments^^ ' 2 ^ 29 ' 34 in 
the linear regime. 

Still, these simplifying approximations oblige us to ne- 
glect a certain number of phenomena, like the interplay 
between nonparabolicity and Coulomb interaction 23 , 
or the scattering toward electron-hole pairs at large 
wavevectors, that is known to be an important fac- 
tor in the thermalization and dynamics of exci- 
ton polariton a 35 ' 36 . Moreover neglecting intrasubband 
Coulomb scattering does not allow us to consider 
polariton-polariton Coulomb interactions, that would ap- 
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pear as an higher order term in the Hamiltonian of Eq. 
(fTT)) . Indeed, the bosonic approximation in Eq. (|Tg|) 
is satisfied only if its expectation value is taken between 
states corresponding to a weak excitation density. Other- 
wise, we will have a correction term in Eq. (|18p . roughly 
proportional to the ratio between the number of excita- 
tions, iV exc , and of electrons, iV e inQw, in such a way that 
this deviation can become important at higher excitation 
densities^ 2 -. 

To go beyond these limitations, we will derive an ef- 
fective quartic bosonic Hamiltonian, taking care of both 
nonbosonicity and Coulomb scattering. In order to do 
so, we will have to calculate the matrix elements of the 
full Hamiltonian in Eq. ([2]) between states with multi- 
ple intersubband excitations. The algebra necessary to 
calculate the coefficients of such an effective Hamilto- 
nian is rather heavy, and we will use the coboson formal- 
ism, originally developed for exciton polariton s 28 ' 37 . The 
reader can refer to the next section for a brief summary 
of such a formalism, highlighting the aspects relevant for 
the algebra of the present paper. 



III. COMPOSITE BOSON COMMUTATOR 
FORMALISM 

Here we generalize the composite boson (coboson) ap- 
proach, which was originally developed in Refs. [13| and 
[28j for the case of excitons. Such a formalism applied 
to intersubband excitations provides a systematic way 
to calculate many-body matrix elements. For simplicity, 
we consider the case where tiqw = 1. The case with 
many quantum wells is more cumbersome, but relatively 
straightforward. Final results are obtained by replacing 
1/JVd by l/(n QW N cl ) in Eq. fE), ([28]) and ©. 

We start by writing Hamiltonian Hp in Eq. @ in the 
appropriate form to use the composite boson commutator 
formalism. In particular it must annihilate the ground 
state. Our main interest will be to calculate matrix ele- 
ments describing scattering events between intersubband 
excitations. We will thus limit ourselves to consider ma- 
trix elements of the form (F\ . . . bi^Hb\, , . . . \F), where 
all intersubband excitation creation operators are on the 
right of the Hamiltonian and annihilation operators on 
the left. It is thus immediate to verify that terms in the 
Hamiltonian containing operators ct k or ci.k in normal 
order, with k outside the Fermi sea, do not contribute to 
such matrix elements, as they annihilate the state. We 
can thus remove them by truncating sums over k, relative 
to the first subband, to the Fermi sea (k < kp). We then 
rewrite Hp in Eq. @ in the following form (omitting 
the ground state energy and non-contributing terms) 

Hp = HKin + -ffRabi + -fflntra + #Dcpob (19) 



where 



H K in = E ^2,fe4,k C 2,k - ^ ^l,feCl,k4 ik , (20) 



Hi n 



1 t>2222 t t 

2 1^ V 1 C 2,k+q C 2,k'-q C 2,k'C 2 ,k 

k,k',q 

+ 9 E ^r mic i,k+q c i.k'-q c i,k' c i, 



H 



k,k',q 

- E n 1221c l,k+ q 4,k<-qC2,k'< k , (21) 
k,k',q 

Depo^— En 12 " K^q 

q 

+ ^0,q 6 0,-q + VqV-q) , (22) 

all sums over the first subband are again truncated to 
the Fermi wavevector and 

nu llk = hw ltk - E ^k-k'i (23) 

|k'|<fc F 

^2,fc = E ^k-k'l- 

|k'|<kf 

In the previous equation we have thus explicitly sepa- 
rated the intersubband and intrasubband terms of the 
Coulomb interaction, putting them respectively into 
-HDepoi and -ffintra- The first term will give the de- 
polarization shifts, while the second one contains the 
electron- hole interaction which lowers the energy of bare 
intersubband excitations- 2 - and is also responsible for the 
scattering between intersubband excitations. Moreover, 
we commuted the electronic operators in such a way 
that now i?Kin and H\ ntra annihilate the ground state 
\F). Writing the Hamiltonian in such a form makes it 
also apparent the renormalization effect of the Coulomb 
interaction over the electronic dispersions, that is the 

same as the one that would be given by an Hartree-Fock 

i 27 so 
approach iisi. 

The starting point of the composite boson formalism is 

the exact commutator between intersubband transitions 

defined in Eq. ^ 



h h 1 



<5m,i^q,q' ^mq t iq ; , (24) 



where D 



mq.zq' 



+ Z)( 2 ) is the deviation from 



bosonicity operator defined by 

n(l) 1 \ " * t 

D = ^2^ m > k ^k+q-q'Cl,k+q-q'Cj jk 



D(2) = ^E^"> k ^k4,k+q' C 2,k+c 



(25) 



To calculate many-body matrix elements involving inter- 
subband excitations, we need to evaluate the following 
commutator 



D 



mq,iq' , u j.q" 



^E{ A ^-<1') (26) 

+ A^(q-q")}4 l ,q' + q»-q' 
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where 



(q) = lTj-^2vm,UVnM+<l l 'ik+< i Vj,k- ( 2 7) 



The term appearing in Eq. (|2l)jl is analogous to an ef- 
fective scattering due to the nonbosonicity of the inter- 
subband excitations (Pauli scattering): two intersubband 
excitations can exchange their electrons to create two ex- 
citations with different wavevectors, the total wavevector 
being conserved. 

With the identities above mentioned, it is possible to 
compute the scalar products of many-excitation states as 



(^1 b n,q'"&m,q&i, q '&j, q » \F) = <J mq ,jq' 5 n q'" , 



' ,01" 



-^^(q-qO + Kq^n.q'")- 

iV e l 



(28) 



From Eqs. (| 13[) and ([22)1 we see that matrix elements of 
-Hrtabi and -ffDepoi can be written as scalar products of 
many-excitation states and can thus be calculated with 
the previous formula. 

In order to evaluate the many-body matrix elements 
of the kinetic part of the Hamiltonian i?Kin, we need the 
commutator 



2,q' 



(29) 



where we have approximated, for simplicity, W2,|k+q| with 
u>2,k- Many body matrix elements of -ffxin can thus be 
calculated by commuting it all the way to the right, leav- 
ing behind only terms in the form of scalar products be- 
tween states with many intersubband excitations. 

To complete the framework, we need to consider the 
matrix elements of -ffintra- We start again by calculating 
its commutator with an intersubband excitation creation 
operator 



H Intra, &j >q 



I>A • l <->- (30) 

3 



where jij is a coupling between bright and dark excita- 
tions 



7i. 



q 



^ 221 A°;?(q), 



(31) 



and Vi t q is called the creation potential operator, 



^ = E [Vq 

Q,m 



V, 



1221 \0,0 



CUQ) (32) 



X ^ni,q+Q 



E C 2,k-Q C 2,k 



+ E (35 1U *B (Q) - Vq 21 ** 

Q , rn 

X6 m,q+QE Cl . k +Q C lk- 



The coefficient 7^ in Eq. (p0|) is a renormalization of the 
intersubband transition energy due to the electron-hole 
interaction. Notice that the creation potential operator, 
V^q, as well as -ffKin and -ff intra, annihilate the ground 
state \F), a property which will be important when cal- 
culating matrix elements. The last commutator we need 
is the one between a creation potential operator and an 
intersubband excitation creation operator, that is 



1 



e 2 n el 



N cl 2e e r 

E£m,i(Q) ,f 

m , n , Q 



(33) 



m,q+Q"n,q'-Q' 



where 



cm = 



\n,j 



t-2222 . \0,0 
1 Q + A m,i 



(Q)A^(Q)/^ 111 

Ko*Zmh 221 > 



is the so-called direct scattering. 

In order to calculate an arbitrary matrix element con- 
cerning -ffintra between many-excitations states, we can 
commute it all the way to the right, leaving behind terms 
that will have the form of a creation potential operator, 
V^q, sandwiched between multiple intersubband excita- 
tion operators. We can then commute creation potential 
operators all the way to the right, leaving at the end only 
many-excitations states scalar products, that we can cal- 
culate using Eq. 



IV. EFFECTIVE BOSONIC HAMILTONIAN 

Our aim is to derive a bosonic Hamiltonian capable of 
describing intersubband polariton-polariton interactions. 
We will thus introduce the bosonic ground state \G) and 
-Bq, the annihilation operator for a bright bosonic inter- 
subband excitation, with in-plane wavevector q satisfying 
the usual bosonic commutation rules 



£q,5q, 



(34) 



The index has been omitted as we are only interested in 
the dynamics of bright excitations. We will then calcu- 
late the coefficients of the bosonic Hamiltonian Hb , im- 
posing that it has the same expectation values that Hp, 
for states with up to two excitations. While the most 
general Hamiltonian, up to n-body interactions, can be 
calculated with the same procedure, here we will limit 
ourselves to consider only resonant two-body interactions 
involving bright intersubband excitations. 

The actual calculation of the coefficients requires some 
rather complex operatorial algebra. For sake of sim- 
plicity, we have thus moved the bulk of the calcula- 
tions in the Appendices, presenting here only the fi- 
nal form of the Hamiltonian and the parameter depen- 
dency of its coefficients. In Sec. IIII1 we have provided 
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the essential ingredients of the composite boson commu- 
tator approac h 22 ' 28 i 37 , that we then use extensively in 
Appendix |E] to calculate the actual coefficients of the 
bosonic Hamiltonian. 



The bosonic Hamiltonian, describing the bright exci- 
tations, and using the quantities defined in Sec. IIIII and 
in Appendix [B] has thus the form 



h b = J2 ( toi2 - ™.°) B l B « + J2 



e 2 n e i 



rl212 

~9 



4eoer.ei 2 (0) q 



2B. 



9=0 



t 5q + (l-C) {b^ + B, 



q -q / 



^ (B q + sL q )(a- q + a q ) 

q 

VEt^^ ^ - *P."»'^)] ^q+P^'-p^'^ + h ' C - 



^QwAoi 



q,q ,p 



1 1 



q,q' p 



(35) 



where the last nonlinear coefficient, describing the scat- 
tering of pairs of intersubband excitations, is 



U 



2 rl212 




(i-C) - 



9=0 



x(0). 



(36) 



In the previous expressions n c \ denotes the electron 
density in each quantum well. A comparison with Eq. 
(fTT)) shows that taking into account the nonbosonicity 
of intersubband excitations and the intrasubband terms 
of the Coulomb interaction resulted in the appearance 
of quartic terms in the Hamiltonian, that makes it pos- 
sible to observe nonlinear polariton-polariton processes. 
In the first and fourth line in Eq. ([33]) and in Eq. (j3"6")l 
coefficients have no dependence over the wavevectors of 
the intersubband excitations. This is justified by the fact 
that these coefficients evolve significantly on a scale of the 
order of the Fermi wavevector kp , which is much larger 
than the photonic wavevectors of interest. There is also 
no spin dependence as one should expect from a naive 
comparison with exciton a 38 i 39 . As explained in Sec. Ill CI 
and Appendix [A] this is due to the fact that all intersub- 
band of interest are made of two fermions with the same 
spin. 

The first quartic term in Eq. (|35|) describes the satu- 
ration of the light-matter interaction. Its physical origin 
lies in the non-perfect bosonicity of the intersubband ex- 
citations and accounts for a renormalization of the Rabi 
frequency at high excitation densities. Its presence is 
easy to understand if one remembers^ that the light- 
matter interaction is proportional to the square root of 
the total electronic population in the first subband minus 
the one in the second subband. As the creation of one 
intersubband excitation corresponds, in average, to the 



promotion of one electron from the first to the second 
subband, we do expect that an increase in the intersub- 
band excitations population will reduce the light-matter 
coupling, consistently with the minus sign in front of such 
a coefficient. 

The second quartic term describes instead the scat- 
tering between pairs of intersubband excitations. It has 
three contributions as seen in Eq. (|36p . The first one 
comes from the intersubband Coulomb interaction. The 
second and third contributions instead are due to the 
intrasubband Coulomb interaction. It is important to 
notice that the two nonlinear coefficients behave differ- 
ently with respect to the number of quantum wells in the 
cavity. While Ml q is proportional to the square root of 
the total number of electrons in the cavity, U depends 
only on the number of electrons in each well. This is due 
the fact that the photonic mode couples to all electrons in 
the cavity while the Coulomb interaction does not couple 
one well with another. 

As we can see, Coulomb interaction and the non- 
perfect bosonicity also give a renormalization effect vis- 
ible in the quadratic part of the Hamiltonian. The 
intersubband energy is renormalized both by Hartree- 
Fock terms and by the Coulomb electron-hole interaction. 
Moreover, the non-resonant terms of the depolarization 
shift are renormalized by a Q < 1 coefficient as a conse- 
quence of the nonbosonicity. Again, the expressions of 
the coefficients can be found in Appendix [B] 

The previous Hamiltonian is adapted to describe op- 
tically pumped experiments, in which the populations of 
the dark states remain negligible. While we will not treat 
this case in the present work, we notice that, in the case 
of electrical pumping, the dark states population can be- 
come significant. Interaction between bright and dark 
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modes is described by the Hamiltonian 



GL 



«#0,q,q' : p 



E 



' n QW N el 

l^O,q,q',p ^ 



where the index i means that the mode is dark. In a mean 
field approach this Hamiltonian can be seen to describe 
energy shifts and saturation effects on the bright modes 
as a function of the dark state population. Coefficients 
in Eq. (1571) can be calculated from a generalization of 
the formula given in Sec. IIIII and in Appendix [Bl 

Passing from the fermionic description Hp to the 
bosonic one Hb, we have kept only quartic terms in 
H B , discarding all higher order terms. This is effec- 



tively equivalent to limit ourselves to the lowest order 
in the excitation fraction iV exc /(nQ W ./V el )2&. We thus ex- 
pect that our theory will correctly describe the polariton- 
polariton interactions in the regime of moderate excita- 
tion densities-^. In order to clarify this point we give here 
an explicit example not involving Coulomb interaction. 
The reader can refer to Appendix [C] for more complex 
and detailed examples. 

We consider a process in which a state with one pho- 
ton and iV exc intersubband exitations is converted into 
a state with no photon and N exc + 1 intersubband ex- 
citations. The matrix element of the Coulomb part of 
the Hamiltonian vanishes in this case, as it conserves the 
number of photons, and we can thus perform the calcu- 
lation exactly, both in the fermionic and in the bosonic 
formalism. We obtain respectively 



J 



(F\a^ T H F bl N - +1 \F) 



\F) (F\ < r + Xq 1^) 



(G\a 1l BZ-H B Bl»-+i\G) 



^/{N cxc + 1)W CXC ! 



m q y/N cxc + 1 \1 - 

m q + 1 ( i - 



N K 



"QW-^el 
A^exc 

A/exc 

2nQ W N ei 



O 



N 2 



(38) 



which are indeed equal to first order in N cxc /(nQwN c \). A. Renormalization of the intersubband transition 



V. NUMERICAL RESULTS 

In this Section we provide the numerical values of the 
relevant coefficients that appear in the bosonic Hamil- 
tonian in Eq. (|55|) . highlighting the dependence of the 
main parameters. 



We start by analyzing the renormalization effect due 
to the Coulomb interaction on the quadratic part of the 
effective Hamiltonian. Performing a Bogoliubov transfor- 
mation on the first line of Eq. (|35|) we obtain the renor- 
malized energy of the intersubband excitation (that we 
should call, more properly, intersubband plasmon3^~— ) 



J 



-ElSBT — 



\ 



7o,o 



e 2 n el 



71212 

_q 



2eoe r ei 2 (0) q 



9=0, 



e 2 n e i 



7-1212 

_1 



2e e r ei 2 (0) q 



(i-0 



(39) 



q=0 



In Fig. ([3]) we plot the dispersion of the intersub- 
band transition energy considering a GaAs quantum 
well of length Lqw = 11 nm (left panel) and Lqw = 
39 nm (right panel) , corresponding to bare transitions of 
140 meV— i and 15 meV— respectively. The solid line de- 
picts the result for the intersubband transition including 



all contributions in Eq. (f3"5j). The dashed line repre- 
sents the result obtained while neglecting the Coulomb 
intrasubband interactions and the nonbosonicity. The 
dash-dotted line is the result obtained as the dashed 
line, but neglecting also the screening of the intersub- 
band Coulomb interaction. Notice that the renormalized 
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FIG. 3: Renormalized energy -Eisbt of the intersubband tran- 
sition for a bare energy hw\2 = 140meV (left panel) and 
15meV (right panel), as a function of the electron density. 
Solid line: result including all the contributions in Eq. (|39[1 . 
namely Coulomb interaction and nonbosonicity. Dashed line: 
result neglecting intrasubband Coulomb interactions and the 
non-bosonicity. Dash-dotted line: same as the dashed line, 
but with no screening of the intersubband Coulomb interac- 
tion. 




02468 10 1 2 3 



n cl (xlO 11 cm~ 2 ) n c i (xlO 11 cm~ 2 ) 

FIG. 4: Thick solid line: effective interaction energy U be- 
tween intersubband excitations including all the contribu- 
tions in Eq. (|36l) f° r = 140meV (left panel) and 
fiwi2 = 15meV (right panel) as a function of the electron 
density. Dashed-line: the first term in Eq. (|36l) correspond- 
ing to the intersubband Coulomb interaction. Thin red line: 
the absolute value of the second term in Eq. (|36|) . namely 
the direct Coulomb interaction. Note that this term is nega- 
tive, thus producing a red-shifted contribution. Dash-dotted 
line: the third term in Eq. (|36p . which is due to the exchange 
Coulomb interaction. 



intersubband energy, -Eisbt, in Eq. converges to the 
bare transition energy hwi2 for vanishing doping. We see 
that the main effect is due to the screening of the inter- 
subband Coulomb interaction. It appears clear that these 
renormalization effects are relevant in the THz regime 
(right panel), but rather moderate for mid- infrared tran- 
sitions (left panel). 



B. Interaction energy between intersubband 
excitations 

Now, we consider the quartic part of the Hamiltonian 
Hb , responsible for the Coulomb scattering of pairs of in- 
tersubband excitations. The interaction energy is given 
by N CKC U/(riQwN e i). In Fig. g] we plot the energy U 
(thick solid line) for a mid-infrared transition (left panel) 
and a THz transition (right panel). The other lines de- 
pict the individual contributions of the three terms in 
Eq. ([55]) (see caption for details). For the considered 
realistic parameters, the interaction energy grows with 
increasing electron doping density almost linearly for the 
mid-infrared case, while it saturates at high densities for 
the case of a THz transitions. Notice that in the case of 
exciton a 28 ' 38 , the direct scattering vanishes in the limit of 
small wavevectors, whereas it is not the case for intersub- 
band excitations when the screened Coulomb interactions 
are considered. This is due the fact that electron-electron 
and hole-hole interactions are screened differently by the 
Fermi sea^. Note that for a doping density in the range 
of a few 10 11 cm -2 , the energy U is of the order of a few 
meV both for the cases of THz and mid-infrared transi- 
tions. The interaction energy can thus reach values close 
to the meV when the density of excitations in the system 
becomes significant. This is rather promising, since as 
shown in the case of exciton-polaritons^ 2 -, very interest- 
ing nonlinear polariton physics occurs when the interac- 
tion energy becomes comparable to the linewidth of the 
polariton modes. For THz polaritons, state-of-the-art 
samples 25 exhibits polariton width in the meV range. 

VI. POLARITON HAMILTONIAN 

In this Section we consider the interactions in the po- 
lariton basis. For simplicity we limit ourselves to the 
first mode of the cavity and we assume that the number 
of electrons is such that the system is in the strong cou- 
pling regime. The diamagnetic interaction in He as well 
as all non-resonant terms in Hb are thus neglected. A 
Bogoliubov transformation of the quadratic parts of He 
and Hb gives the expression of the polaritonic operators 

\PiJaJ \aq ft, J \a^J 

where pu q and pl q are polaritonic operators of the upper 
(UP) and lower branch (LP) respectively, while a q and 
ft, are the Hopfield coefficients. Note that results can be 
generalized to the ultrastrong coupling regime (vacuum 
Rabi frequency comparable to intersubband transitions) 
by considering the Bogoliubov transformation involving 
the antiresonant light-matter interaction terms^. Here, 
for simplicity, we consider the case of moderate values 
of the vacuum Rabi frequency. Because of the quartic 
terms in Hamiltonian Hb polaritons interact with each 
other through their matter part and can scatter. In the 
following, we will focus on the lower branch. Replacing 
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operators a q and B q by their expression, the two-body 
interaction between LP polaritons is found to be 

#LP-LP = \ E „ q ' q '^ PLq+pPLq'-pPLq'PLq>( 41 ) 



2 , "-QW^cl 



where 



Vq.q'.p = Mq'-p^'K^^q^)- (42) 



As an example, we now consider the case where the 
lower branch is pumped at a wavevector q p so that the 

system is in the state Pl q 8XC | G) / \J iV cxc ! . This Hamilto- 
nian allows us to describe single-mode (Kerr) and mul- 
timode (parametric) coherent non-linearities^. For the 
parametric case, one has to consider interaction channel 



qp 



t^ex 

-p p Lq P 



\G), (43) 



where pairs of polaritons scatter from the pumped mode 
into signal-idler pairs. As for the case of exciton- 
polaritons we expect that the maximum efficiency of this 
parametric processes is achieved when the energy conser- 
vation condition is fulfilled^. A mean-field approach of 
the proble m 42 ' 44 shows that the matrix element between 
the initial and the final states is the relevant quantity to 
consider and has to be compared with the lifetime of the 
excitations. For high pump intensity i.e. N cxc ^> 1 this 
matrix element is 



M, 



-V a 



q p p — n at qp<qp<p - 



(44) 



As discussed in Sec. IV Bl and shown in Fig. SI for THz 
polaritons nonlinear interaction energies of the order of a 
few meV (thus comparable to THz polariton linewidths) 
can be achieved, thus paving the way to a very interest- 
ing coherent nonlinear physics for this kind of composite 
excitations. 



VII. CONCLUSIONS 

In this paper, we have presented comprehensively a mi- 
croscopic theory for the manybody physics of intersub- 



band cavity polaritons. Using a composite boson com- 
mutator approach, we have derived the manybody inter- 
action matrix elements involving (screened) Coulomb in- 
teraction processes and Pauli saturation effects. We have 
also derived an effective bosonic Hamiltonian reproduc- 
ing the matrix elements in the two-excitation manyfold 
calculated by the composite boson approach. As a first 
application, we have calculated the renormalization of 
the intersubband transition frequency at the mean-field 
level and determined the strength of the interactions be- 
tween intersubband excitations. We have also derived 
the renormalization of the light-matter coupling due to 
the Pauli saturation effects and considered the interac- 
tion physics in the polariton basis. Our results show that 
significant intersubband polariton-polariton interactions 
occur, especially for transitions in the THz regime. Our 
work paves the way to promising future studies of nonlin- 
ear quantum optics with intersubband cavity polaritons. 
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Appendix A: Spin and quantum well number 



In this appendix we give the explicit notations for Eq. 
([3]) and Eq. (|12[) with quantum well and spin indexes. 
Electronic wavefunctions are localized in quantum wells 
and we neglect electronic tunneling from one well to an- 
other. We thus neglect Coulomb interaction between 
electrons in different wells, which are sufficiently apart. 
Eq. Q and (|12l) can then be written with all indexes 



J 



«QW 



-ffRabi 
HqouI 



0)t JJ) 



j=l k,o-,/j 

»QW 

E E ^(4iU4i CT +4iU4i.)(«q+« t - q ) 

j=l k,q,(T 

1 nqw 

iV V t/a^V M t r U) t M jj) 

2 / i / j q jU,k+q,<r °i/,k'— q,o-' °i/',k',o-' L /i',k,cr' 



(Al) 
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and 



wave vectors) 



y/n QW N e i 



1QW 

k,<T 



„co t 



y 0,j,k L 2,k+q,cr c l,k 



„,(A2) 



where i^oj.k = 0(feF — for all j, k is a wavevector such 
that k < k F , a, a' € {!,{} and ft, (j,' ,v,v' G {1,2}. Bright 
intersubband excitations are thus linear superposition of 
pairs of fermions with the same spin. We could generalize 
Eq. (| A2|) by allowing the two fermions to have different 
spins but the resulting collective excitation would be dark 
and thus not relevant if we consider only polariton. 



Appendix B: Coefficient calculation 

In order to determine the coefficients of the effective 
bosonic Hamiltonian in Eq. (|35p . we impose that it has 
the same matrix elements as the fermionic one in Eq. @ 
over the Hilbert space spanned by one and two excitation 
states. That is we impose that 



{g\t^h b tI„tI„,\g) 



(F\t^P^H F t\„t\,„ \F) 



\/^q,q'^q" 



(Bl) 



q, -Bq], tq £ 



i,6o,q], P± is the projec- 



where T q e [I, < 
tor over the subspace orthogonal to the ket^i and iV q)q / 
and ^4q, q ' are the normalizations of the orthogonalized 
states. 

The constraints in Eq. (|B1[) . together with the fact 
that we decided to consider, beyond quadratic terms, 
only the quartic ones, are enough to uniquely define Hb- 
The calculation of the coefficients of Hb thus reduces to 
the calculation of matrix elements of the fermionic Hamil- 
tonian H F , that we can perform exploiting the commu- 
tator approach outlined in Sec. IIIII 

We thus start from a general bosonic Hamiltonian 



q q 

+ ^fin,(B q + Bl q )(o_ q + ot) 
q 

q,q ,p 



■-> ^ n^Ni B ™ B «-* B « B * 

q,q :P 



(B2) 



We also introduce the following five matrix elements 
(each one implicitly dependent on all the relevant 



d x = (F\b ,qH F bl^\F) 



(F\ b .qH F a\ \F) 
(F\ HFblJ^ \F) 

(F\ 6 ,q + p60,q'-p^&J,q^q \F) 

(F\ 6o, q+P &o,q'-pP±tfi4 )q 4, q \F) ■ (B3) 



d- 2 
ds 
di 



The relations between the coefficients of Hb in Eq. (|B2|) 
and the matrix elements in Eq. (|B3|) can be found with 
some algebra to be 

K = di 
2Q 



y/1 - 2/nQwN e l 

d 2 



(<5 P ,o + V-q,p) ^« 



G 



q-q ,q -q-p 



'N, 



q+p,q'-p 



.4 



q+p,q'-p 



(B4) 



and 



2 (<5 p ,o + <*q'-q,p) K ' 



2U 



ffiq.q'-^q+p.q'— p 
Aj : q'^q+p,q'-p 



d 5 , (B5) 



where, for consistency, terms under the square root have 
also to be developed to the first order in l/(riQvjN e \). 

We have thus reduced the determination of the coeffi- 
cients of Hb to the calculation of the five matrix elements 
in Eq. (|B3|) . Without loss of generality, in the following 
calculations, we will assume that the excitations wavevec- 
tors are small compared to the Fermi wavevector, which 
implies that the Pauli scattering involving two bright ex- 
citations, defined in Eq. (|27j) . can be approximated by a 
Kroneker delta 



(P) 



(B6) 



1. Coefficient di 

The only contributing terms are the kinetic and 
Coulomb part of the Hamiltonian Hp. The matrix el- 
ement is 



12 — 7o,o + 



2 7-1212 
e n cl lq 



2eoe r ei 2 (0) q 



(B7) 



9=0 



12 



2. Coefficient di 



to the normalization. 



Only the light-matter part of Hp contributes and we 
obtain 



3. Coefficient d% 



(B8) 



4. Coefficient d4 

As for the coefficient d,2, only the light-matter part of 
Hp contributes, giving 



This term comes from the intersubband Coulomb in- 
teraction and is given by 



e 2 n e i 



7-1212 



2e e r ei2(0) 



(1-0, (B9) 



9=0 



where 



rl212 
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c = -E- 



9 =0 



Q 



' r A°;°(Q). (BIO) 



The coefficient Q is then 0^(1 + 1/tiqw-Wci)- As the num- 
ber of electrons is large we can neglect the correction due 



J 



Ml q I 5 p o + ^p,q'-q 



"QW^cl 



(Bll) 



The last term, at the origin of the first quartic term in 
Eq. (1551) . is due to the nonbosonicity. 



5. Coefficient d$ 

The expression of the coefficient is given by the follow- 
ing relations (to first order in l/(<T.Qw-^ei)) 



d, = ((F\ & , q +p&0, q <-p - (F\ &0,q+p&0,q'-p&S,q'&0,q \ F ) <*1 &0,q&0,q') Hp^ ^ \F) 



(B12) 



(F\ bo^+pbo.q' -pHpbl ^bl ^ \F) = 2 (h£u 12 - 70,0) (^ P ,o + <5 P , q '- q - 

2 7-1212 / 

Op.O + Op.q' 



£0er£l2 



9=0 



4-2C 



I 



/n e ] 



2eoe r 



^(0) , 



(B13) 



where 

3f(p) 



^ n iV cl Q £ (Q) °< m 

m , n , Q 



A°'l(Q-p)C°o(-Q)(B14) 



In Eq. (|B13|) the long wavelength limit has been taken 
because the variation of the coefficients is significant only 
on a scale of the order of the Fermi wavevector kp. 



Appendix C: Transition rates and lifetimes 

In this Appendix we calculate the transitions rates 
for scattering of pairs of intersubband excitations us- 
ing both the effective bosonic Hamiltonian approach and 
the fermionic formalism developed by M. Combescot and 



coworkers^, showing that the two approaches give ex- 
actly the same results 



1. Fermionic case 

In this paragraph we use the composite boson commu- 
tator formalism to calculate the transition rate between 
an initial state \ipi) and a final state |^/), as well as the 
lifetime of the former, using a first-order, time-dependent 
perturbation theory. The calculation is equivalent to a 
Fermi golden rule as explained in Ref. (37j . 

The particular event we want to describe is the scat- 
tering of an initial pump beam of arbitrary intensity into 
a signal and an idler mode. We will thus consider initial 
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and final states to be proportional respectively to the kets As in Appendix [B] Pi is the projector over the subspace 

orthogonal to \ipi) and F t verifies 



l^> * bl N - \F) 



and 



IV>/> oc bl^bl^-pbl^ 2 \F) 



The transition rate can then be written as^ 



(CI) 



(C2) 



\F t {E)\ 2 = ^S t (E), 



(C5) 



(C3) 



where where St converges to the Dirac delta function for long 

times. The orthogonal projection in Eq. (|C4|) can be 
= F t (H F - {^Hf^PxHf^) . (C4) calculated as 



P^C"" \F) 



-N cxc (hio 12 - 7o ,o) E^^tqC"" 1 l F > 



A^exc(-^e: 



1) 



ATexc(iV, 



1) e 2 n cl y-, £"o°(Q) fc t 



2n QW A cl 2e e r ^ Qe(Q) n ' q_CJ °' q 



2n QW A^ el 



A CXC (A CXC -1) e 2 n G i 
2riQ\vA c i eo£ r £i2 



^A™ °(Q-q)6t iQ 6L, 2q _ Q C-^ 2 |P) 



m , Q 



r 



(C6) 



The previous expression is a first order in the pertuba- Ft, i.e we do as if intersubband excitations were exact 
tions. We can thus consider only the zero-th order in eigenstates of Hp, obtaining 



J 



(F\ 6^- 2 6 0jq _p6 , q+p F t (iJf - (VM H F !</>;)) « F t (0) (F| Vq-P&o.q+P- 



(C7) 



Taking into account the normalization we thus obtain the to be of the form 

transition rate 1 1 ^-^ 



T 



2tt N exc (N e 



1 



n 2 N 2 

"•QW iv ol 



-U 2 5 t (0) + O 



"QW^cl 



(CIO) 



.(C8) 

where it is implicitly assumed that the summation is re- 
stricted to small wavevectors. To obtain this result we 
made F t act at zero-th order and we used the fact that 
With the same notations, the contribution of the bright and dark intersubband excitations are not reso- 
many-body physics to the lifetime of the initial state is21 nant. The counterintuitive 1/2 factor comes from the 



5; = (i/t 



(C9) 



The second term can be calculated using the same 
method as for the transition rate and is found to con- 
tribute only to higher orders in the perturbation. After 
some rather cumbersome algebra, the first term is found 



overcompleteness of the composite boson basis. 



2. Bosonic case 

In this paragraph we calculate the same quantities as 
in the previous one but we assume that intersubband ex- 
citations are bosons and that their dynamics is described 
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by Hamiltonian Hb- We can therefore use a traditional 
Fermi golden rule. 

The initial and final states in this case are 



m ex Bl N - \G) 



(Cll) 



and 



cx Bl +p Bl_ p Bl N -- 2 \G). 
The transition rate is thus given by the formula 

2n N exc (N exc - 1) 2 



h n 2 N 2 



U 2 6 t (0), (C12) 



which is the same as in Eq. (|C8j) up to third order in 
N exc /(nQwN e i). As all bright intersubband excitations 



are resonant the lifetime of the initial state is the sum of 
the transition rates 



(C13) 



where the summation is again restricted to small 
wavevectors. A comparison with Eq. (|C10|) shows that 
this method underestimates the true lifetime by a fac- 
tor two. This is coherent with the results in Ref. [37$ . 
that show how an effective Hamiltonian approach that 
correctly calculates transition rates needs to take into 
account an ad hoc factor 1/2 when calculating lifetimes, 
due to the overcompleteness of the composite boson ba- 
sis. This is of course simply implies a renormalization of 
the composite boson density of states. 
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